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In heavy ion collision experiments as well as in neutron stars, both baryon and isospin chemical potentials are 
different from zero. In particular, the regime of small isospin chemical potential is phenomenologically important. 
Using a random matrix model, we find that the phase diagram at non-zero temperature and baryon chemical 
potential is greatly altered by an arbitrarily small isospin chemical potential: There are two first order phase 
transitions at low temperature, two critical endpoints, and two crossovers at high temperature. As a consequence, 
in the region of the phase diagram explored by RHIC experiments, there are two crossovers that separate the 
hadronic phase from the quark-gluon plasma phase at high temperature. 



1. INTRODUCTION 

Heavy ion collision experiments are important 
for our understanding of the strong interaction 
at nonzero temperature and density. In heavy 
ion collision experiments, both the baryon and 
the isospin densities are different from zero. The 
time between the formation of the fireball and its 
freezeout is so short that only the strong interac- 
tions play a significant role and baryon number as 
well as isospin are conserved. Therefore, it is phe- 
nomenologically worthwhile to study the effects of 
a nonzero isospin chemical potential, fj,j, on the 
QCD phase diagram at nonzero temperature, T, 
and baryon chemical potential , [Ib- However, 
most studies at high temperature and nonzero 
density have been restricted to cases where either 
fii or [Ib is zero. 

For jiB 7^ and /zj = 0, the phase diagram is 
very rich. At low T and high fj,s, the ground state 
is believed to be a color superconductor pQ. If T is 
increased and fiB is decreased, a first order phase 
transition separates the hadronic phase from the 
quark gluon plasma phase. If /is is further de- 
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creased, this first order critical line stops in a crit- 
ical endpoint. At lower there is a crossover 
between the hadronic phase and the quark gluon 
plasma phase. These results are based on differ- 
ent models |2I3| . as well as on exploratory lattice 
studies at low chemical potential 0]. 

In the case of (j.b — and /xj ^ 0, the fermion 
determinant is real and traditional lattice meth- 
ods can be applied. In this case, the phase dia- 
gram is very rich as well. At low T, an increase in 
[ii above the pion mass leads to a supcrfluid phase 
with a pion condensate. At low T, the phase tran- 
sition between the hadronic phase and the pion 
condensation phase is second order and has mean 
field critical exponents. If T is increased, this 
second order phase transition becomes first or- 
der. Therefore there is a tricritical point in the 
phase diagram. These results were found using 
both lattice simulations |S] and effective theories 
UJ. At high T and low /if, a crossover separates 
the hadronic phase from the quark gluon plasma 
phase. There also might be a critical endpoint 
and a first order phase transition at high T when 
Hi is increased 0. 

We use a Random Matrix model as a schematic 
model for QCD to study the phase diagram at 
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nonzero T, fig, and pi [7|. This model has been 
previously used to study QCD at Pbj^O and pj = 
[3]. We then analyze possible consequences for 
heavy ion collision experiments. 

2. RANDOM MATRIX MODEL 

Random matrix models were introduced in 
QCD to describe the correlations of low eigenval- 
ues of the Dirac operator [S]. It was shown that 
for large matrices these models are equivalent to 
the mass term of the chiral Lagrangian that is 
uniquely determined by the symmetry of QCD 
[H]. Therefore, in the chiral limit, Random Ma- 
trix Theory provides an exact analytical descrip- 
tion of the low-lying Dirac spectrum. The idea is 
to replace the matrix elements of the Dirac oper- 
ator by Gaussian random variables subject only 
to the global symmetries of the QCD partition 
function. The temperature and chemical poten- 
tials enter through external fields structured ac- 
cording to these symmetries. 

For two quark flavors of mass m\ and m 2 , the 
partition function of our model reads 
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where Wi = W + fl + m, and Wi = + ft 1 * - m 
are n x n matrices, with f2 = diag(iT, —iT). The 
elements of the random matrix W are complex. 
First, we express the determinant in Q as a 
Grassmann integral. Then, we perform the Gaus- 
sian integration over W. Third, the resulting 
four-fcrmion interaction is decoupled by means 
of a Hubbard-Stratonovich transformation at the 
expense of introducing mesonic degrees of free- 
dom. Finally we perform the Grassmann inte- 
gration, and the partition function can be 
mapped onto |7j 



Z = JvA exp(-£(A,A t )), 
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where A is an arbitrary complex 2x2 matrix, and 

C = nG 2 Tr(A - M^)(A^ - M) - | Tr log Q+Q, (3) 

with M = diag(mi, m^), and 

A iT + [1q + pih 
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with I 3 = diag(l, —1). This is an exact mapping. 

In the large-n limit, this partition function can 
be evaluated by a saddle point approximation. 
To solve the saddle point equations, we make an 
Ansatz for the matrix A based on the symmetry 
of the partition function. Since we have two inde- 
pendent chemical potentials, one for each quark 
flavor, the chiral condensates are not necessarily 
equal. Furthermore, at high enough pi, we ex- 
pect a pion condensate. We thus assume that 
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In this parameterization, the chiral condensates 
are given by (uu) — G 2 (ai — mi), (Jd) = G 2 (cr 2 — 
7712), and the pion condensate by ^((u^d) — 
(J75U)) = G 2 p. We thus get an effective poten- 
tial that can be studied in the usual way. 

We limit ourselves to the case mi = m 2 = m. 
We are particularly interested in the phase dia- 
gram in the p,B-T plane at pti small enough so 
that the superfluid phase is never reached, be- 
cause it corresponds to the conditions of heavy ion 
collision experiments. The phase diagram in the 
/^s-T-plane for zero p,i has been studied in pjj. In 
the chiral limit, the chiral restoration transition 
extends as a second order line from the ps = 
axis, changes order at a tricritical point, and in- 
tersects the T = axis as a line of first order tran- 
sition. For nonzero quark mass, the first order 
transition ends in a critical point, and the second 
order transition becomes a crossover. Figure ^ 
shows the phase diagram in the /is-T-plane at 
finite quark mass mG — 0.1 for zero isospin chem- 
ical potential, pi = 0, and for pjG — 0.1. We 
observe that the first order curve splits into two 
first order curves that are separated by 2pjG. 
This can be understood as follows. Below the 
threshold for pion condensation, the free energy 
separates into a sum over the two flavors. For 
pi = 0, the chiral phase transition lines for both 
flavors coincide. A nonzero pj breaks the flavor 
symmetry, and the first order transition lines for 
the two flavors split and shift in opposite direc- 
tions. The temperature of the critical endpoints 
is not affected by pi. 
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Figure 1. Phase diagram in the (J-b-T '-plane for 
quark mass roG = 0.1 and [ii as shown. First or- 
der transitions are denoted by full curves, and 
crossovers by dotted curves. The condensates 
a i = (uu) and a 2 — (dd) are omitted where <C 1. 

3. CONCLUSIONS 

We have used a Random Matrix model for 
QCD at nonzero T, /i_b, and fii. We have found 
that in the region of high T and small fin , an arbi- 
trarily small [ij greatly alters the phase diagram 
in the hb~T plane: There are two crossovers, two 
critical endpoints, and two first order phase tran- 
sition lines that separate the hadronic phase from 
the quark gluon plasma phase ? . This could 
have important consequences for heavy ion colli- 
sion experiments, since they are done at \ib =/= 
and [ii 7^ 0. If our Random Matrix model gives 



an accurate description of the phase diagram, it 
might also be interesting to use different isotopes 
in heavy ion collision experiments in order to vary 
Hi at constant hb- These results have been con- 
firmed by other models, albeit with some con- 
straints in one case ^U]. Finally, recent lattice 
studies at small nonzero /zb can also be used at 
nonzero [Ib cmd fij, and will provide an important 
test for our results. 
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